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We present a parameter- free method of extraction of the electron magnetic circular dichroism 
spectra from energy- filtered diffraction patterns measured on crystalline specimen. The method is 
based on multivariate curve resolution technique. The main advantage of the proposed method is 
that it allows extraction of the magnetic signal regardless of the symmetry and orientation of the 
crystal. In other words, our new method overcomes difficulties caused by complexity of dynamical 
diffraction effects. 



I. INTRODUCTION 

Electron magnetic circular dichroism^ (EMCD) is the 
name for a transmission electron microscopy analogue of 
the x-ray magnetic circular dichroism (XMCD) 25 . Both 
techniques allow studying magnetic properties of sam- 
ples via an excitation of a core electron into unoccupied 
valence states. Due to sharp atom-specific energy lev- 
els of the core electrons, this makes both EMCD and 
XMCD element-selective. Transmission electron micro- 
scope (TEM) brings an additional advantage of a supe- 
rior lateral resolution, which reaches in principle sub- 
Angstrom range. 

Despite an intense initial activity in EMCD, which 
led to improvements of spatial resolutiorP^, advances 
in theoretical descript ion^*^, benchmarking quantitative 
measurement J 13 * 15 " 17 ! and first application s 18 " 2 1 \ the ex- 
periments are still primarily in the method development 
stage. The main reasons are high demands on quality of 
sample preparation, low signal to noise ratio (SNR) and 
complications due to dynamical diffraction effects. 

While the first problem can only be approached from 
the experimental side, the other two can be aided from 
the theory side. So far the best control over signal to 
noise ratio one can obtain by measuring energy-filtered 
TEM datacubes, using the telefocus mode^. This way 
one has a two-dimensional set of spectra, which can be 
processed independently one-by-one to provide a rich 
statistics for error control. It seems to be generally ac- 
cepted that the most suitable experimental geometry is 
the 3-beam orientation, in which one can use the infor- 
mation from all four quadrants in the diffractio n plane by 
means of the double-difference (DD) procedur d 16 | 17 | 22 [ 

Here we present a method, how to overcome the re- 
quirement of the 3-beam orientation and allow to ex- 
tract the EMCD signal from arbitrary crystal orientation. 
The method is based on advanced statistical techniques, 
namely the multivariate curve resolution (MCR^P. In- 
troduction of this method into the EMCD field means 
that the requirement of a particular crystal orientation 
is not anymore required - thus essentially removing the 
obstacles caused by dynamical diffraction effects. 



In the following text we will describe the principles and 
assumptions of the MCR method (Sec. |n|, method of 
construction of theoretical datacubes used in this article 



(Sec. Ill) and the results obtained (Sec. IV). 



II. MULTIVARIATE CURVE RESOLUTION 



The present multivariate analysis is based on 
the MCR-alternating least-square (ALS) algorithm^, 
whereby we adopted a faster, accurate and more robust 
algorithm, called the non-negative matrix factorization 
(NMFpHISl The main improvement of the present NMF 
is to include smoothness constraints, which are often en- 
forced to regularize the computed solutions in the pres- 
ence of noise in the data. 

The single EEL spectrum extracted from the datacube 
at a position corresponding to a single pixel is supposed 
to be a linear combination of two kinds of spectral com- 
ponents, each consisting of nonmagnetic and either clock- 
wise or counter clockwise chiral components with dif- 
ferent weights, depending on the position in reciprocal 
space. The objective of MCR-ALS is to extract pure 
component spectra with as few assumptions about the 
data as possible by repeatedly applying the least-square 
fit until the residual is converged within the experimental 
noise level. 



As reported previousl y 27 * 28 1 , the MCR technique is be- 
set with two inherent difficulties: one is that the number 
of components underlying cannot be determined a pri- 
ori; the other is that the result obtained by means of 
the MCR technique is not unique. The first problem 
disappears in the present EMCD data. The second one 
is usually compensated for by including auxiliary con- 
straints on the solutions. We found that the additional 
smoothness constraint dramatically improved our second 
problem. The MCR solutions were computed using the 
NMF method as described in^. 
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III. CONSTRUCTION OF DATACUBES 

We have tested the MCR technique on theoretical dat- 
acubes. The advantage of this approach is that we know 
precisely the spectral shape of the EMCD signal and its 
distribution throughout the diffraction plane. 

As a benchmark system we chose a bcc iron crystal. 
The electronic structure was calculated in density func- 
tional theory with local density approximation using the 
code WIEN2k 2 ^ with over 100 basis functions per atom 
and 40000 k-points in the 1st Brillouin zone, which en- 
sures well converged bandstructure. Based on that, we 
calculated the MDFFs 30 using the A = 1 approximation^. 
The dynamical diffraction effects were treated using the 
Bloch's waves formalism described in the same article. 
We used the automatic Bloch-wave summation technique 
MATS 31 with criterion 10 -5 , which leads to highly con- 
verged calculations in a computationally efficient way. 

We have calculated a datacubes in two orientations: 
1) A slightly mis-oriented 3-beam orientation so that the 
extraction of the EMCD signal would be non-trivialp^ 
yet favorable for the DD procedure. The chosen two ex- 
cited Bragg spots were G = ±(200) and Laue circle cen- 
ter (-0.095, 0.078)G. 2) A two-beam orientation with 
an excited Bragg spot G = (110) and Laue circle center 
(t^0). Here we are rather far from the symmetrical 3- 
beam orientation and we can expect difficulties in EMCD 
extraction due to asymmetry 32 . 

In both cases, the energy-filtered diffraction patterns 
span an area from — 1.2G to 1.2G along the systematic 
row and from — 0.7G to 0.7G in the perpendicular di- 
rection, with step 0.025G in both directions. The en- 
ergy range for calculation of MDFFs was from zero up to 
15 eV above the Fermi level with step 0.1 eV, separately 
for both L 3 and L 2 edges. That means a starting dataset 
of dimensions 97 x 57 x 151 per edge and thickness. Sub- 
sequently, the Ls and edge spectra were first shifted 
in energy up to 708 eV and 721 eV, respectively, which 
corresponds to experimental values of onset of the L 2 ^~ 
edges in iron. The spectra were then smeared by 0.7 eV 
and 1.2 eV wide Lorentzian for L3 and L2 edge, respec- 
tively, to account for the finite lifetime of the excitations. 
Finally, they were overlapped into a single datacube con- 
taining both edges. The resulting datacube has an energy 
range from 695 eV to 745 eV with a step of 0.05 eV. Sam- 
ple thicknesses were chosen as 10 nm, 20 nm, 30 nm and 
40 nm. 

On top of the £2,3 signals, we added a double-step 
background signal modeled as a sum of two step functions 
positioned at 708 eV and 721 eV, respectively, broadened 
by the same Lorentzians as the edge signal, i.e., 0.7 eV 
and 1.2 eV, respectively. The relative ratio of heights of 
the two step functions was set to 2:1. The relative inten- 
sity of the post-edge background signal was fixed to 1/3 
of the intensity of the broadened nonmagnetic part of the 
L3 edge peak. These choices reflect the following assump- 
tions: 1) the dynamical diffraction effects have the same 
influence on the double-step background signal intensity 



as on the 2p — ^ 3d transitions, 2) the double-step back- 
ground signal originates from transitions to non-polarized 
states, and, 3) the life-time of the transitions in question 
is the same as for £2,3 edges. 

Small deviations from these assumptions are possible 
in experimental situations, for example in 33 the polar- 
ization of the s states is discussed, with its influence on 
the extracted orbital to spin moment ratio from XMCD. 
However, our main aim is to test an accuracy and robust- 
ness of the MCR, therefore a study of influence of such 
effects goes beyond the scope of present manuscript. 

Additionally, for each of the four thicknesses, we have 
also generated datacubes with random noise. Two dif- 
ferent noise levels were generated. In particular, we have 
multiplied the datacubes by a constant factor, so that the 
L3 peak intensity J max becomes 1000 or 10000 counts, 
respectively, and all the other values were scaled accord- 
ingly. Afterwards, at all p ixels and en ergies a random 
noise value from range ±y/I(q x , q yj E) was added. This 
allowed us to compare the performance of MCR with DD 
procedure at various noise levels. 



IV. RESULTS 

Here we summarize and discuss our results. First we 
show the raw magnetic signal spectrum and its distri- 
bution in the near-3-beam orientation-as a reference for 
testing the quality of DD and MCR extraction. The 
datacubes are then analyzed by DD procedure and in 
the next subsection by MCR. In the last subsection we 
present an analysis of the 2-beam orientation datasets in 
a more concise form. 



A. Net EMCD signal 




FIG. 1. Distribution of the L3-edge integrals of EMCD sig- 
nal of bcc iron crystal in the diffraction plane at 4 different 
thicknesses. The energy integral of the L3 edge signal at the 
transmitted beam is normalized to one. The area marked in 
the 1st quadrant by a dashed line was used in DD method for 
extraction of EMCD signal. 
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FIG. 2. Averaged EMCD spectra (full lines) and their energy 
integrals (dashed lines) of bcc iron crystal in the diffraction 
plane at 4 different thicknesses. The energy-integrated EMCD 
spectrum is scaled so that its value p at 718 eV equals to 1. 
Post-edge value q enters the sum-rule expressions, see text for 
details. 

As indicated in the previous section, for a tilted 3-beam 
orientation there is no symmetry operation, that would 
allow an exact extraction of the EMCD signal from the 
datacube. The DD procedure was devised to minimize 
the errors introduced due to misorientation and it was 
successfully applied ir i 16 * 17 i Here we work with a theo- 
retical datacube and thus we have access to the EMCD 
signal by construction. Namely, we calculated the initial 
dataset in two steps, first only considering the real parts 
of the MDFF and second, considering only the imagi- 
nary parts of the MDFF. It is the latter that contains 
a net EMCD signal, neglecting the tiny magnetic dipole 
tempo] gy a ppiyi n g the Lorentzian broadening and over- 
lapping L 3 and L 2 edges, as discussed above, we obtain 
datacubes containing net magnetic signal. This serves as 
a reference dataset, against which we can test the results 
of the DD and MCR analysis. 

Fig. [I] shows the distribution of the EMCD signal in 
the diffraction plane. In the figure we show an energy 
integral of the EMCD at L3 edge and the color bar range 
is relative to the transmitted beam intensity at L3 edge. 
We can see that the overal intensity of the EMCD signal 
decreases with thickness and the pattern of its distribu- 
tion becomes more complicated due to stronger dynam- 
ical diffraction effects at larger thicknesses. The maxi- 
mum intensity of EMCD signal reaches around 0.3% of 
the transmitted beam intensity - the major reason of 
difficulties with obtaining high SNR in experiments. 

Fig. [2] summarizes net EMCD spectra averaged over 
the diffraction plane. The averaged spectra were ob- 
tained by summing all EMCD spectra, choosing the sign 
of L3 magnetic contribution as negative. The plot con- 
firms that the magnetic signal decreases with thickness, 
as anticipated from Fig. [I] By scaling the individual 



EMCD curves to the same peak magnitudes, the shapes 
of spectra become practically identical. This is because 
the magnetic properties in our model are based on elec- 
tronic structure of bulk iron regardless of sample thick- 
ness considered in the dynamical diffraction calculation. 

In quantitative EMCD experiments^ 1111 , the E MCD 
spectrum is analyzed by sum rule expression d 12 * 13 * in or- 
der to extract the ratio of orbital to spin magnetic mo- 
ment. We will apply these sum rules as a sensitive test 
of the accuracy of EMCD spectrum extraction. The sum 
rules state 

m x = 2 f L Av(E)dE + f L Av(E)dE _ 2q 

m s 3 J L Aa(E)dE - 2 J L Aa(E)dE 9p - 6q U 

where q is an energy integral of EMCD over both edges 
and p is an energy integral over the L3 edge only. With- 
out loss of generality we can rescale the EMCD spectrum 
or its integral so that p = 1. Then the mi/m s ratio is a 
function of q only, allowing easy visual comparison of the 
energy integrals in the post-edge region. 

In this manuscript, the integral over L3 edge will mean 
an energy range from 695 eV to 718 eV. There is a cer- 
tain degree of ambiguity in this choice, however we have 
checked that none of our conclusions do substantially de- 
pend on this choice. 

Performing the energy integral of the scaled (p = 1) net 
EMCD spectra leads to a value of q = 0.079 which gives 
mi/m s = 0.019. This is in good agreement with the DFT 
values rrij = 0.048/i# and m s = 2.26/i#, i.e., rrij /m s = 
0.021. The error can be attributed to energy dependence 
of the radial wavefunctions in matrix elements 34 . 



B. Double difference procedure 

In Ref. 16 we have introduced the DD procedure on 
post-edge normalized datacube. This method efficiently 
cancels the systematic errors introduced by slight mis- 
orientations from a perfect 3-beam orientation^. 

When applying DD, we first normalize the datacube 
in the post-edge region, so that the energy integral over 
the last 2.5 eV is equal for all spectra in the datacube. 
(We will discuss later the dependence of results on cho- 
sen normalization range.) Then for every energy- loss, a 
mirror image of the lower diffraction half-plane is sub- 
tracted from the upper diffraction half-plane (1st differ- 
ence) and the lower half-plane is removed from further 
considerations. Finally, we subtract a mirror image of 
the left quadrant from the right quadrant (2nd difference) 
and we remove the left quadrant from further considera- 
tions. I.e., the dataset is reduced by DD to one quarter 
of its original size. Now a spectrum at every pixel cor- 
responds to an approximation of the EMCD spectrum, 
assuming an efficient cancellation of non-magnetic com- 
ponents by the DD. An averaged EMCD spectrum can 
be constructed either from the whole quadrant or some 
selected region, as was done irP^. 
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We will apply this method of extraction of EMCD 
spectra to our theoretical datacubes (both noise-less and 
noisy ones) to a) test the accuracy of the DD procedure 
on realistic model data with known magnetic component 
and b) to compare its performance to the new method 
of extraction of EMCD based on MCR discussed in the 
next subsection. 

First we have extracted averaged EMCD spectra from 
the whole quadrant. At selected four sample thicknesses, 
from the noise-less datacube we obtained stable mi/m s 
value of 0.029. For the noisy datacube with peak- value 
I max = 10 4 the span of mi/m s values was 0.023-0.028 
and for the most noisy datacube (7 max = 10 3 ) the span 
was from —0.013 to 0.052, which is a large error. 



stable results for the noise-less data: we obtained q = 
0.121 ±0.003, which gives mi/m s = 0.029. Note that this 
value over-estimates the expected value 0.019 based on 
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FIG. 3. EMCD spectra and their energy integrals for bcc 
iron crystal at 4 different thicknesses obtained by the double- 
difference procedure from a region marked in Fig. [I] The 
energy- integrated EMCD spectrum is scaled so that its value 
p at 718 eV equals to 1. The inset zooms in on the post-edge 
region of the energy integrals, q. The values of q enter the 
sum-rule expressions for evaluation of mi/m s ratio, see text 
for details. Full, dashed and dotted lines correspond noise- 



less, I u 



10 and 10 , respectively 



While the result is quite satisfactory for noise-less or 
low-noise (7 max = 10 4 ) datacubes, for the most noisy dat- 
acube the accuracy of mi/m s extraction is quite poor. 
The reason is that in DD the non-magnetic part of the 
spectra is almost removed and in areas of zero or low 
magnetic signal we are adding basically only noise and 
eventual anisotropy signal. That was also the reason for 
using a limited integration window around an area of the 
strongest magnetic signal instead of the whole diffrac- 
tion plane irPH That introduces some additional param- 
eters, such as position, shape and size of the integra- 
tion window. Here we have used a rectangular region 
around the so-called Thales circle positiorP with dimen- 
sions 0.5G x 0.45G, similarly as i n 1 see Fig. [l] We 
have tested that small variations in the position and size 
of the region do not qualitatively change the results. 

Use of the limited integration window led again to 



raw EMCD evaluation in Sec. IV A Our analysis shows 
that this offset is a systematic error and that it comes 
from the post-edge normalization procedure itself, see 
Appendix [B] The error is approximately 0.01 for our 
dataset and is almost independent of the value of mi / m s 
ratio itself. The magnitude of the error depends on the 
chosen broadening model and range of energy-integrals 
over edges, particularly in the post-edge region. 

With increasing noise in datacubes the accuracy of DD 
relatively quickly drops down, see Fig. [3j For the dat- 
acube with peak-value 10 4 the DD procedure provided 
values mi/m s = 0.029-0.035. At highest noise level - 
datacube with peak-value 1000 - the DD procedure pro- 
vided large spread of values mi/m s — 0.007-0.044 (note 
the order of magnitude spread), yet somewhat improved 
in comparison to the evaluation based on the whole quad- 
rant. 



C. MCR method 

Results obtained by MCR method are described here. 
First we discuss the efficiency of decomposition of the 
spectra in terms of maps of coefficients and the residu- 
als and then we summarize the extracted mi/m s ratios. 
Note that in MCR method we do not use any integration 
window as in DD procedure. The reason is that MCR effi- 
ciently uses all data available to strengthen the model of 
both magnetic and non-magnetic spectral components, 
i.e., areas with low EMCD signal help to improve the 
model for non-magnetic spectral component, which in 
turn effectively improves the accuracy of extraction of 
the magnetic one in areas with stronger EMCD signal. 
In this sense, use of a limited integration window can be 
counter-productive for MCR and tests performed on our 
datasets have confirmed that (not shown). 

The MCR procedure outlined in Sec. [TT] produced sta- 
ble fits for all 12 tested datacubes. Every such fit re- 
sults in isolation of two spectral components along with 
maps of corresponding positive coefficients throughout 
the diffraction plane. The distribution of the EMCD sig- 
nal is then obtained as a difference of coefficients of the 
two spectral components. These maps are summarized 
in rows 2, 4 and 6 of Fig. [Z] for noise-free, J max = 10 4 
and J max = 10 3 datacubes, respectively. These maps can 
be compared to theoretical relative dichroic maps shown 
in the top row and an excellent agreement is found. The 
only differences are due to noise, but the position and 
shape of regions of positive or negative magnetic signal 
is recovered with high accuracy. 

A more quantitative measure of the quality of the fit 
is expressed in terms of maps of residui. The value 
of a residuum at a particular pixel (q x ,q y ) is defined 
as v^J^fe,^,^)] 2 , where AI(q x ,q y ,Ei) represents 
the difference between the actual and fitted value of the 
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FIG. 4. Maps of distributions of EMCD signal and residui 
of MCR fits. Row 1 shows distribution of the raw magnetic 
signal in post-edge normalized datacube (cf. Fig. [I]). Row 
2 shows map of coefficients of the EMCD signal extracted 
by MCR on noise- free datacube. Row 3 shows an associated 
map of residui. Rows 4-5, and 6-7 show the same as rows 
2-3, but for the noisy datacube with peak level 10 4 and 10 3 , 
respectively. Each column corresponds to different thickness, 
as indicated in the top row. 



pixel at energy E { . For the noise- free datacube the map 
of residui shows negligible values demonstrating that the 
model fitted by MCR represents data with high accu- 
racy. Small deviations are found only at positions of 
Bragg spots, but even there the error is negligible (see 
the color scale). 

For the noisy datacubes, residues necessarily can't be 
zero. Overall, SNR is approximately 3.2 times higher for 
the datacube with i ma x = 10 4 compared to the datacube 
with J max = 10 3 . That can be explained by the nature 
of Poissonian noise: the datacube with J max = 10 4 has 
10 times stronger signal, the noise increases by factor 
of a/10 and in consequence SNR improves by factor of 
10/a/IO ~ 3.2 - in agreement with the scales of the maps 
of residui. 

The shape of the noise distribution can be explained 
via the normalization. Since the intensity of the signal 
in the corners is weaker than around the transmitted 
beam or Bragg spots, SNR is lower in the corners. Af- 
ter the post-edge normalization of spectra, the noise in 
low-signal regions becomes amplified. 

It is interesting to note, that noise is the dominant 
source of residui in our noisy datacubes - making the fit 
more accurate around peak positions, while it is less ac- 
curate in low signal areas. For the noise-free datacube it 
is the other way around, see row 3 in Fig. [4] and compare 
with rows 5 and 7. It indicates that the spectral shape 
is slightly different at Bragg spots, compared to regions 



in between. Since the Bragg spots form only a fraction 
of the area of diffraction plane, in presence of zero noise 
they show the highest residual signal. It would be inter- 
esting to understand the nature of this very weak, but 
systematic residual signal, however that goes beyond the 
scope of the present study. 
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FIG. 5. EMCD spectra and their energy integrals of bcc 
iron crystal at 4 different thicknesses obtained by the MCR 
method. See caption of Fig. [3] and text for details. 

Fig. [5] summarizes the quantitative analysis of the 
EMCD spectra extracted by MCR method. For the 
noise-less datacube we obtain in agreement with DD pro- 
cedure a value of q = 0.120 with a spread of less than 
0.001, leading to a value mi/m s = 0.029 - again contain- 
ing the systematic error originating from the post-edge 
normalization, Appendix [Bj For the noisy datacube with 
/ max = 10 4 we obtained mi/m s ratio ranging from 0.019 
to 0.036. This spread is comparable to the DD procedure, 
however, contrary to DD procedure, in MCR the range 
of values is centered around mi/m s = 0.029 extracted 
from noise-less datacube. Finally, for the noisy datacube 
with J max = 10 3 we obtained mi/m s values ranging from 
0.017 to 0.051. 

This shows that in a near-3-beam orientation, which is 
optimal for application of the DD procedure, the MCR 
method is at least as efficient in extraction of EMCD 
spectra as the DD procedure. In the next subsection 
we will compare performance of both methods far from 
3-beam orientation, where DD is less applicable. 

In the Table U we summarize and extend the results of 
both DD and MCR. We have included an error bar, which 
comes from an analysis of the obtained mi/m s ratio as 
a function of post-edge normalization interval, namely 
from 1 eV up to 5 eV. In general, we have observed the 
following behavior: a) for too narrow intervals, partic- 
ularly in noisy datasets, the extracted mi/m s ratio can 
vary strongly-therefore we picked 1 eV as the smallest 
range; b) for too large intervals the obtained mi/m s ra- 
tio slowly and monotonously grows, because the fraction 
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TABLE I. Summary of mi/m s extraction in near-three-beam 
orientation. Net values are obtained directly from averaged 
EMCD spectra. Then DD and MCR methods were applied on 
noise- free and noisy datacubes. The level of noise is expressed 
via the maximum intensity J max at the L3 peak. 
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see Eq. (B2). 



magnetic signal within a quadrant, see Fig. [6] By fine- 
tuning the size and shape of the virtual aperture it might 
be possible to improve the results of DD procedure. On 
the other hand, this shows the potential caveats of this 
method. Note that in RefP^ only an up-down differ- 
ence was applied to extract the EMCD signal in the two- 
beam orientation. Since we are far from three-beam ori- 
entation, one might wonder whether a simple difference 
would be more accurate. We have tested this, however 
the results with simple up-down difference were worse or 
comparable with DD procedure (not shown). 

In the same way as in the near-3-beam orientation, 
the MCR method demonstrated its robustness also here. 
The EMCD spectra were easily extracted and the maps 
of magnetic signal well corresponded with the expected 
calculated distributions, see Fig. [6] The mi/m s ratios 
were extracted with high accuracy, except for the noisiest 
datacubes at thicknesses 20 nm and 30 nm, see Table \H\ 
However, even there the MCR method provided values 
closer to expected mi/m s = 0.028 than DD procedure. 



CONCLUSIONS 



D. Two-beam orientation 
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FIG. 6. Same as Fig.^] but in two-beam orientation with G = 
(110). The integration window tested with DD procedure is 
indicated in the top row. 

Raw EMCD spectra in the two-beam orientation are 
expectedly very similar to those in Fig. [2j providing also 
the same mi/m s ratio. The DD procedure did not per- 
form very well, especially for 20 nm and 30 nm datasets, 
see Table [TTJ The reason of poor results at these thick- 
nesses can be understood due to low relative strength of 
the magnetic signal and also due to varying sign of the 



We have developed a new method for extraction of 
the EMCD spectra using multivariate curve resolution 
(MCR). We have shown that with this method we can 
extract the EMCD spectrum with the same or higher 
precision than with the double-difference procedure. The 
main advantage of MCR method is its generality - it 
does not require presence of any symmetry planes. The 
MCR method is therefore applicable in any geometry, 
not only in the near-3-beam orientation required by 
double-difference procedure. In principle, by use of MCR 
method we avoid the obstacles caused by dynamical elec- 
tron diffraction. Our results will simplify future analyses 
of the experimental EFTEM datacubes and extraction of 
the EMCD spectra. 



TABLE II. Summary of mi/m s extraction in two-beam orien- 
tation. Net values are obtained directly from averaged EMCD 
spectra. Then DD and MCR methods were applied on noise- 
free and noisy datacubes. The level of noise is expressed via 
the maximum intensity J max at the L3 peak. 
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30 


nm 


0.0181 


10 4 


0.03 ±0.19 


0.032 ± 0.002 








10 3 


-0.29 ±0.10 


0.070 ±0.019 








00 


0.029 ±0.002 


0.029 ±0.001 


40 


nm 


0.0180 


10 4 


0.023 ±0.005 


0.022 ± 0.004 








10 3 


0.032 ±0.014 


0.025 ± 0.009 
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Appendix A: Brief description of MCR-ALS method 

The datacube is rearranged into a two dimensional ma- 
trix form, each column corresponding to a single spec- 
trum extracted from a single pixel of the datacube. Given 
mxn matrix X representing the datacube, mxk matrix 
of concentration profiles C, and nxk matrix of pure com- 
ponent spectra S, where k is the number of components, 
the following relation holds 



X = CS T +E, 



(Al) 



where the superscript j denotes the transpose of the ma- 
trix. E is the residual matrix with the data variance 
or statistical noise unexplained by CS T . The goal is to 
obtain physically meaningful C and S. The data set con- 
sists of a set of intensities J(x, E) at the spatial position 
(x, y) and energy-loss E. The data can thus be expressed 
by the three dimensional array X(m £C , m y , n), whose vari- 
ables respectively correspond to the pixel (channel) coor- 
dination numbers of the x and y directions and energy- 
loss E. However, the spatial coordinates can be equiv- 
alently treated and thus the data is essentially bilinear, 
expressed by the array X(m,n), where m = m x x m y . 
Now each row in the matrix S T on the right hand side of 
Eq. ( Al ) stands for a pure component EEL spectrum and 



if each row vector of S T is normalized in its intensity (and 
the total core-loss cross sections are not very different for 
the component spectra), the ra-th value of the k-th col- 
umn of the matrix C denotes the relative composition of 
the k-th component spectrum (k-th row of S T ) at the ru- 
th spatial coordination. This is a two-way bilinear model 
and the alternating least-square (ALS) method is usu- 
ally utilized as a solution of a straightforward regression 
problem where initial estimates for C and S T are gen- 
erated initially using either random numbers, eigenvalue 
decomposition, or dissimilarity criterion. 



Appendix B: Influence of normalization on 
quantitative results 

Schematically we can decompose a spectrum into non- 
magnetic part TV, magnetic part M and double-step back- 
ground signal B: 

a(E, q x , q v ) = N(E) + M(E) + B(E) 



In the post-edge (PE) region the following relations 
can be extracted from our datacubes: 



Bp E = Nl 3 , max/3 

M PE « -0.002M Ls , max 
N PE w O.OOliV. 

.max 



(Bl) 
(B2) 
(B3) 



These values naturally depend on the chosen broadening 
model and the definition of the PE area. 

Generally the magnetic signal is at most 10-15% of the 
non-magnetic signal at L3, i.e., 



\M L 

q ,max I < 0.157V L , 



(B4) 



therefore the post-edge intensity of the spectrum in our 
model datasets is 



cr(q x ,q y )p E ~0.001N L3 



Nl 3 



0.002M 



Z/3,max 

(B5) 

where we used relation Bpe = iVL 3)max /3 applied in con- 
struction of our datacubes. 

To simplify the notation, we'll introduce the following 



v(qx,q y )pE = C(i±a) 

C = N L „ 0.001 
a = 0.002|M max |/C 



(B6) 
(B7) 
(B8) 



Using this notation and Eq. (B4), the value of 



a < 0.002 x 0.15A/L 3 , max /C « 0.0001 <C 1 
and therefore 

[C7(l±a)]- 1 «(l T a)/C7. 



(B9) 



That will be used below in normalization of spectra. 

Now we will turn our attention to energy integrated 
spectra. We need separate energy-integrals over the L3 
and L2 edges, respectively. We take a difference of two 
PE- normalized spectral integrals to obtain edge- 

integrals of the approximate magnetic signal Aa: 



Aa 



(Nl 3 

2 
C 



[Ml 3 



M Ls + B L 
M L3 + B L 

1 a(N Ls -i 



■N L2 -M L2 +B L2 )(l + a)/C 
■N L2 +M L2 +B L2 )(l-a)/C 



B L3 )-M L2 +a(N L2 +B L2 )} 



Here, M^ 2 3 ,^ 23 and B^ 23 stand for edge-integrals of 
magnetic, non-magnetic and background signals, respec- 
tively. Note that in the difference of the PE-normalized 
spectra there are terms originating from non-magnetic 
and background signal. These can introduce a system- 
atic error into mi/m s evaluation. 



8 



Let's try to estimate a magnitude of this error accord- 
ing to sum rule, eq. 0: 

m _ 2 [Ml 3 + a(N L3 + B L3 )} + [-M L2 + a(N L2 + B h2 
m s ~ 3 [M L3 - 
2 Mr, 



- B L3 )] - 2[-M L2 + a(N L2 + 5 L2 
a(7V Ls +7V L2 +5 L3 +5 L2 ) 



3M L 

777^ 

m. q 
2 
3 

7774 



-a- 



net 



■ 2M L2 + a(AT L 
3 



,3) 



Si 



r L 3 + 2M L2 
M L2 )(B Ls -2B L2 ) 



1 — a 



+ 2M i2 )2 
-Bl 3 - 2B L2 



Mr 



2Mr 



2a N L2t3 



3 Mr 



2M L 



where we assumed Nl 3 — 2Ni 2 <C Bl 3 — 2Bl 2 . 

The size of the last term can be estimated using 
Eqns. fB8) and pi|: 



2a N L2t3 



Bl 2 



0.004|M„ 



N + B 



3 M L3 +2M L2 N L3)max (l+ 0.003) 3M 



(BIO) 



where we assumed that M L:i ~ M^ 2 = M. Assuming 
that the shapes of the magnetic peak and non-magnetic 
one do not vary much throughout the diffraction plane 
we get M oc M max and iV, B oc A^ 3j max , i.e., the ratio 
^max/^L 3 ,max x (N+B)/M is a constant and its value for 
our dataset is approximately 2.24, which gives an error 
estimation 



0.004 x 2.24 « 0.01 



(Bll) 



Therefore this term alone explains the difference between 
expected 0.019 and obtained 0.028 (MCR and DD) values 
of mi/m s ratio. 

The second error term is proportional to mi / m s . Ratio 



2B t 



0.002|M max | 2B L2 



M L3 + 2M L2 



Nl 3 



M 



(B12) 



is again approximately constant throughout the diffrac- 
tion plane, because M ex M max and Bj J2 3 oc NL 3 ^ max . Its 
value is around 0.002 x 5. 7 = 0. 012, which after multipli- 
cation by mi/m s gives a negligible correction 0.00023. 
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